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1. Introduction 

Remarkable progress achieved in the description of phase transitions and critical 
phenomena due to application of the renormalization group pQ (RG) ideas leads 
sometimes to the conclusion that all principal work in the field has been done, especially 
if universal features of criticality are addressed. This is certainly not true if one goes 
beyond a description of the para-to-ferromagnetic phase transition in a standard N- 
vector model j2], which belongs to the universality class of the 0(iV)-symmetric 4 
theory. To give an example, a study of realistic systems often calls for account of a 
single-ion anisotropy j3] , structural disorder |Hj or frustrations [I] which might essentially 
change the critical behaviour. An account of this new physics still remains a challenging 
problem. 

The theoretical RG description of the above mentioned systems requires field 
theories of complicated symmetry with several couplings. In theory, the critical point 
corresponds to the fixed point of the RG transformation. An accessibility of the fixed 
point, along with the (non-universal) initial conditions of the RG flow is governed by the 
so-called marginal field dimensions. These are universal and together with the critical 
exponents and amplitude ratios constitute intrinsic features of criticality. It is well 
established by now that the universal features of a 2nd order phase transition in the 
3d TV-vector model are not sensitive to the single-ion cubic anisotropy if iV < iV™ b 
jsJj. They are neither changed upon a weak quenched dilution by a non- magnetic 
component if only 0[ N > N^ 11 . Currently, there exists a good agreement between 
the numerical values of the marginal dimensions iV£ ub , calculated in numerous RG 
and MC studies j3J EJ HO]- However, this is not the case for the frustrated systems, 
where even more marginal dimensions have been found and their numerical values are 
still under discussion jH [TT] . 

The problem we want to raise in this report concerns the critical behaviour of the 3d 
frustrated spin systems with noncollinear order. The most common physical realization 
of such systems are stacked triangular antiferromagnets (the examples are given by 
CsMnBr 3 , CsNiCl 3 , CsMnI 3 , CsCuCl 3 , VC1 2 , VBr 2 ) and helimagnets (Ho, Dy, Tb, 
(3 — M11O2). In the former case, the noncollinear order is caused by the frustrations due 
to the triangular geometry of the underlying lattice, whereas in the latter one it is due 
to the competition of ferro- and antiferromagnetic interactions. Currently, there exists 
a large literature devoted to the subject, which results from more than twenty years 
long studies 4]. However, neither experimentally nor theoretically has an unanimous 
conclusion been drawn so far about the nature of the phase transition into the ordered 
state in these systems. Important physical quantities which are under discussion are 
the marginal dimensions of the models. In particular, when the model is generalized 
to describe iV-component vectors ^2] (physical systems mentioned above correspond to 
N = 2, 3), one finds a marginal dimension N C3 below which the phase transition is of first- 
order, whereas for N > N C:j it is of second-order. Several variants of the perturbative 
RG expansions [SJ HB1 UH US1 UE| and various truncations for the Wilson-like non 
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perturbative RG (NPRG) equations [T71 [TBI HH] give different numerical estimates 
for N C3 (see Table UJ). However, they all agree, that such a dimension (along with two 
marginal dimensions more, N C2 and N Cl , see below) exists. In our study, we aim at 
performing a thorough analysis of these marginal dimensions by means of the pseudo-£ 
expansion [20]: the technique which is known to provide the most accurate results in 
the 3d perturbative RG approach. 

The rest of the paper is organized as follows: in the next section we formulate 
the model we are interested in and obtain the expansions for its marginal dimensions, 
section 01 is devoted to the numerical estimates on their basis, section Ogives conclusions 
and outlook. 

Table 1. Marginal dimension N C3 obtained within different RG methods. See the 
text for details of the methods. 



fixed 3<i, resummed 
3 loops 6 loops 
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2. The model and the pseudo-e expansion 

An effective Hamiltonian of the model of frustrated magnets with an iV-component 
order parameter reads [T2] : 

n = J &x{\[p<kY + (V0 2 ) 2 + mM + &)} + + ^ + 

In (HJ), mo,UO)^o are Dare mass and couplings and 4>i = 4>i{ x ) are ^-component vector 
fields, representing the cosine and sine modes associated with the spin ordering. The 
noncollinear (chiral) ordering occurs for u, v > 0. 

We sketch the fixed points (FPs) picture retrieved in the previous RG studies [T2] . 
For the high field dimensions iV > N C3 four FPs exist: the Gaussian, u — v — 0, G, the 
Heisenberg u ^ 0,v = 0, 0(2N) symmetric, H (both unstable for the space dimension 
d — 3), and two nontrivial FPs u / 0,n ^ (with u, v > 0), chiral and antichiral C + 
and C-. The FP C + is stable and governs the chiral 2nd order phase transition. With a 
decrease of N, the FPs C + and C_ merge at iV = N C3 and disappear: only unstable FPs 
G and H are present for iV just below N C3 . Note that Pelissetto et al. and Calabrese et 
al. 0IE] have claimed that, once resummed, the six loop /5-functions obtained directly 
in d = 3 exhibit a new root — a new fixed point — below an extra marginal dimension 
of the field estimated at iV ~ 5.7. According to these authors, this fixed point is neither 
analytically related to the Gaussian fixed point in d = 4 nor to the fixed point found at 
large iV in d = 3. It should therefore be non perturbative, although it is found within 
the perturbative framework. Thus, they claim that below this new marginal dimension, 
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and in particular for N = 2, 3, the transition is again of second order. Note also that 
this scenario disagrees with the results obtained within the NPRG approach for which 
the transition is (very weakly) of first order for iV < N C3 . We will return to this point 
at the end of our report. As iV is further decreased, the nontrivial FPs C + and C_, 
existing above N C3 reappear at iV < N C2 , but now in the u > 0, v < quadrant and thus 
do not describe the chiral phase (they both have complex coordinates for N between 
N C3 and N C2 ). Finally at iV = N Cl one of the nontrivial FPs merges with the FP H 
and, with further decrease of N, passes to the quadrant u > 0,v > 0, still remaining 
unstable. The above picture is supported both by the perturbative RG approach (e- 
expansion accompanied by subsequent resummation or by a conjecture about the 
series behaviour resummed expansion in terms of renormalized couplings in a 3d RG 
scheme 00112], 1/N expansion [HUE!) and the non-perturbative RG ff\ \T7 \ IT8 | IT9"]. 

Discrepancies between the values of the marginal dimensions N Ci obtained so far 
within the perturbative RG (see Table [T] for N C3 ) to a great extent are because the series 
used for their analysis were rather short fHfl. IT4*| EE Ej and in general are known to be 
asymptotic at best pQ. So it is very desirable to perform an estimate of the marginal 
dimensions on the basis of the expansions which, on the one hand would be of the 
highest order and, on the other hand, would possess better convergence properties. As 
we show below, these two goals are reached by applying the pseudo-e expansion to the 
six- loop d = 3 RG /5-functions obtained for the effective Hamiltonian (JJJ in the Ref. 

The method consists in introducing an auxiliary parameter (r) into the /3-functions 
which allows to separate contributions to the FPs from the loop integrals of different 
order [20] • This is achieved by multiplying a zero- loop term by r and obtaining FP 
coordinates and, subsequently, all FP quantities of the theory as series in r with final 
substitution r = 1. Starting from the six- loop d = 3 RG /3-functions of Ref. jS] we get 
the expansions for the fixed point coordinates and derive for the marginal dimensions: 

iV C 3 = 21.797959 - 15.620635 r + 0.262060 r 2 - 0.150930 r 3 - 

0.039165 r 4 - 0.029721 r 5 , (2) 

N C2 = 2.202041 - 0.379365 r + 0.202166 r 2 - 0.084951 r 3 + 

0.092744 r 4 - 0.097961 r 5 , (3) 

N C1 = 2.0 - 0.666667r + 0.145212 r 2 - 0.094836 r 3 + 

0.099757 r 4 - 0.112325 r 5 . (4) 

Expansions (J2])-(J1|), derived directly at d = 3, can be compared with the 
corresponding e = 4 — ci-expansions jH] : 

N C3 = 21.80 - 23.43e + 7.088e 2 , (5) 
N C2 = 2.202 - 0.569^ + 0.989e 2 , (6) 
N C1 = 2 - e + 1.294e 2 . (7) 

Formulas P ]) -({4* |) take into account three orders of perturbation theory more than the 
highest available e = 4 — <i-expansions ©-(JH). Moreover, comparing and 
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one sees that the expansion coefficients in the pseudo-e: series decay much faster and 
one may expect to get more convergent results on their basis. And indeed this is the 
case as we will see in the next section. 



3. Numerical estimates of the marginal dimensions 

The field theoretic RG expansions are known to have zero radius of convergence and 
different resummation techniques are used to make numerical estimates on their basis 
P]. Here, we make use of the Pade-analysis [21] to make an analytic continuation of the 
expansions for r = 1. On the one hand already this simple technique allows us to show 
essential features of the pseudo-e expansion behaviour, on the other hand it allows to 
determine numerical values of the marginal dimensions with a sufficient accuracy. The 
results for the pseudo-e expansion series (J2J)-(jI!) are given below in the form of Pade- 
tables. There, a result of an [M/N] Pade approximant is represented as an element of 
a matrix with usual notation, e.g. the first row gives results of the mere summation of 
the series: 



N, 



C3 



N, 



(-'2 



N, 
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1.751 


7.220 








2.202 


1.823 
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1.984 
1.962 


1.966 
1.977 
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0.385 
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1.985 


2.012 
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1.333 
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1.453 


1.421 


1.432 


1.458 
1.421 


1.032 
1.105 

1.436 


1.431 
1.431 


1.431 


1.446 


1.432 






1.415 









6.126 
1.318 



(8) 



(9) 



(10) 
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Results shown by fractions indicate, that the Pade approximant contained a pole for r 
close to 1 (the denominator shows the value of r for which the pole is obtained). These 
tables can be compared with the analogous tables for the e-expansions ©-((Zj): 

21.80 -1.630 5.458 
X. = 10.507 3.812 . (11! 

7.505 

2.622 



Nr. 



Nr 



2.202 
1.750 
2.514 

2 

1.333 
1.813 



1.633 
1.994 

1 

1.564 



2.294 



(12) 



(13) 



One certainly sees that the convergence properties of the pseudo-£ expansion are better 
in comparison with the e-expansion (cf. the convergence of the results along the main 
diagonal and those parallel to it: there the Pade analysis is known to provide the most 
reliable data |21j). One more feature of the expansions for N c . is evident when one 
compares tables IjlOjl : whereas the central elements of the table give a firm 
estimate for N C1 : [2/2] = [3/2] = [2/3] = 1.431, such a stable behaviour is not found in 
the corresponding Pade tables ©, © for N C3 , N C2 . Obviously, this different behaviour 
is connected with the different origin of the marginal dimensions N Cl from the one side, 
and N C3 , N C2 from the other side. Indeed, the dimension N C1 corresponds to merging of 
the non-trivial and Heisenberg FPs after which the non-trivial FP continuously passes 
to the other quadrant of the u — v plane whereas dimensions N C3 , N C2 correspond to 
the coalesce and disappearance of two non-trivial FPs (see discussion at the beginning 
of section EJ) . 

To make the numerical estimates on the base of Pade-tables (|S ]) -([10 p we proceed 
as follows. For N C3 we take on the main diagonal the highest Pade approximant 
with an estimate [2/2] = 6.23 and suppose that the deviations from an account of 
higher-order terms will not exceed the difference [2/2] — [1/1] = 0.21. For N C2 we 
take the highest obtained estimates [3/2] = [2/3] = 1.99 considering a confidence 
interval as [3/2] — [2/2] = 0.04. Subsequently, for N Cl the central value is given by 
[3/2] = [2/3] = [2/2] = 1.43 with a confidence interval [2/2] - [1/1] = 0.02. Finally, we 
get for the marginal dimensions: 



N C3 = 6.23(21), N C2 = 1.99(4), N Cl = 1.43(2). 



(14) 



The above estimates include within the error bars all elements of corresponding Pade- 
tables except of the inverse approximants [0/N] (and the approximant [5/0] of 0) and 
therefore the confidence intervals in (|14j) are rather overestimated. 

Comparison of our estimate for N C3 with the perturbative RG data of Table ^ 
supports recent estimates E] N C3 ~ 6. We also suggest that essential difference 
between this estimate and the numbers obtained within e- and 1/iV-expansions 
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|14[ IT51 Ibo] is because the last have not been estimated with comparative accuracy 
which was caused in particular by shortness of corresponding series [22] ■ Available so 
far estimates of N C2 are due to the resummation of three-loop massive RG expansions 
[T3] and of the e 2 expansion jHj N C2 = 1.96 and N C2 = 2.03(1), correspondingly 
[22] • Together with the symmetry arguments ^3] providing N C2 > 2 our estimate 
suggests that the value of N C2 should be located very close to 2. In particular this 
means that corresponding scenario of appearance of the pair of non-trivial FPs which is 
governed by this marginal dimensions might not be found in numerical calculations 
for N = 2. Dimension N C1 has its counterpart ^31 as the marginal dimension of 
the N- vector model with a single-ion cubic anisotropy: iV Cl = N^ nb /2, see section 
[0 The last has been estimated by different methods. In particular the Pade-Borel 
resummation of the series (|1J) gives [2j the number coherent with the other data jH] 
N Cl = N° uh /2 = 1.431(3). Comparison of this number with our estimate ()14j) based on 
much less elaborated technique supports the reliability of chosen here scheme. 

4. Conclusions 

The numerical values of the marginal dimensions that we obtain represent clearly 
an improvement of the preceding determinations performed both by the perturbative 
methods and the NPRG one. Once again the pseudo-e expansion turns out to be 
very accurate and constitutes probably a new way to analyze the critical behaviour of 
3d frustrated magnets. However, the principal question about the order of the phase 
transition in these systems for N = 2, 3 still remains open. Obviously, our studies 
are in coherence with the FP picture of the NPRG and perturbative RG approaches, 
where no stable FP are found in the region N C3 > N > N C2 (however the difference 
between numerical value of N C3 and typical numerical values found in the NPRG studies 
[3 HZl UH UH| calls for a more detailed analysis). Nevertheless, the existence of the FP 
found recently in Refs. Ej and claimed to describe criticality of these systems can 
neither be supported nor rejected by our perturbative approach. We think that one of 
possible ways to shed light on this problem is to try to follow an evolution of this FP 
with change of d in order to understand its origin at the upper critical dimension d = 4. 
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Note added. After the completion of this work [23], we learned about e 4 -results for 
N C3 , N C2 and a six-loop pseudo-e result for N C3 [21] . The five-loop e-expansion improves the 
e 2 data for N C3 = 6.1(6) but leads to the unphysical conclusion N C2 = 1.968(1) < 2. The 
pseudo-e expansion of Ref. |24| for N C3 coincide with our formula (J2J), but we report higher 
numerical accuracy. 



On the criticality of frustrated spin systems with noncollinear order 



8 



References 

[1] Zinn- Justin J 1996 Quantum Field Theory and Critical Phenomena (Oxford: Oxford University 
Press); Kleinert H and Schultc-Frohlindc V 2001 Critical Properties of <fi - Theories (Singapore: 
World Scientific) 

[2] See e.g. Holovatch Yu (ed) 2004 Order, Disorder, and Criticality (Singapore: World Scientific) to 
appear 

[3] For a recent review of the critical behaviour of systems with complicated symmetry see: Pelissetto A 
and Vicari E 2002 Phys. Rep. 368 549. Disordered systems are further addressed in: Folk R, 
Holovatch Yu and Yavors'kii T 2003 Physics- Uspekhi 46 169 [Uspekhi Fizicheskikh Nauk 173 
175] 

[4] For reviews, see: Diep H (ed) 1994 Magnetic Systems with Competing Interactions (Singapore: 
World Scientific); Delamotte B, Mouhanna D and Tissier M 2003 Nonperturbative 
renormalization group approach to frustrated magnets Preprint cond-mat/0309101 and Refs. 
El E| IB] for more recent theoretical and MC studies 
[5] Pelissetto A, Rossi P and Vicari E 2001 Phys. Rev. B 63 140414(R) 
[6] Calabrese P, Parruccini P and Sokolov A I 2003 Phys. Rev. B 68 094415 
[7] Tissier M, Delamotte B and Mouhanna D 2003 Phys. Rev. B 67 134422 
[8] Itakura M 2003 J. Phys. Soc. Jap. 72 74 

[9] Folk R, Holovatch Yu and Yavors'kii T 2000 Phys. Rev. B 62 12195; Erratum: Phys. Rev. B 63, 
189901(E) (2001) 

[10] Bervillier C 1986 Phys. Rev. B 34 8141; Dudka M, Holovatch Yu and Yavors'kii T 2001 J. Phys. 

Stud. 5 233; and unpublished 
[11] Parruccini P 2003 Phys. Rev. B 68 104415 

[12] See: Kawamura H 1988 Phys. Rev. B 38 4916 and references therein 

[13] Antonenko S A and Sokolov A I 1994 Phys. Rev. B 49 15901 

[14] Antonenko S A, Sokolov A I and Varnashev K B 1995 Phys. Lett. A 208 161 

[15] Pelissetto A, Rossi P and Vicari E 2001 Nucl. Phys. B 607 [FS] 605 

[16] Gracey J A 2002 Nucl. Phys. B 644 [FS] 433 

[17] Zumbach G 1993 Phys. Rev. Lett. 71 2421 

[18] Tissier M, Delamotte B and Mouhanna D 2000 Phys. Rev. Lett. 84 5208 

[19] Tissier M, Delamotte B and Mouhanna D 2001 Int. J. Mod. Phys. A 16 2131 

[20] The pseudo-£ expansion technique was introduced by B.G. Nickel, (unpublished), see Ref. [19] in 

Le Guillou J C and Zinn-Justin J 1980 Phys. Rev. B 21 3976. For an application of the pseudo-e 

expansion in determination of marginal dimensions see Refs. [!|J I10| 
[21] Baker G A, Jn and Graves-Morris P 1981 Pade Approximants (Reading, MA: Addison- Wesley) 
[22] The confidence interval ±0.01 found in the three-loop studies of Ref. ^2] should be considered 

rather as an accuracy estimate of the calculation scheme, since the six-loop studies within the 

same perturbation technique ® shift e.g. the central value for N C3 from 3.91 to 6. 
[23] Delamotte B, Holovatch Yu and Ivaneyko D 2003 On the critical behaviour of stacked frustrated 

triangular antiferromagnets 2 nd Int. Conf. "Physics of Liquid Matter: Modern Problems" 12-15 

September, 2003 Kyiv, Ukraine" ed L Bulavin (Kyiv) p 72 
[24] Calabrese P, Parruccini P 2003 Five- loop e expansion for 0(n) x 0(m) spin models Preprint 

cond-mat/0308037 v2 



